International Journal of Power Electronics and Drive System (IJPEDS) 
Vol. 10, No. 4, December 2019, pp. 1914~1922 
ISSN: 2088-8694, DOI: 10.11591/ijpeds.v10.14.1914-1922 O 1914 


A model for predicting photovoltaic module performances 


Nadia Bouaziz!, Arezki Benfdila?, Ahcene Lakhlef® 


' Departement of Electrical Engineering, University of Tizi-Ouzou, Algeria 


2,3 Departement of Automatics, University of Tizi-Ouzou, Algeria 


1,2,3 Micro and Nanoelectronics Research Group, Faculty of Electrical Engineering and Computer Sciences, 


University Mouloud Mammeri, Algeria 








Article Info 


ABSTRACT 





Article history: 


Received Mar 4, 2019 
Revised Apr 14, 2019 
Accepted Jun 24, 2019 


Keywords: 


Conversion efficiency 
Electrical parameters 
I-V characteristics 


The present paper deals with the development of a simulation model for 
predicting the performances of a solar photovoltaic (PV) system operating 
under current meteorological conditions at the site location. The proposed 
model is based on the cell equivalent circuit including a photocurrent source, 
a diode, a series and shunt resistances. Mathematical expressions developed 
for modeling the PV generator performances are based on current-voltage 
characteristic of the considered modules. The developed model allows the 
prediction of PV cell (module) behavior under different physical and 
environmental parameters. The model can be extended to extract physical 
parameters for a given solar PV module as a function of temperature and 
solar irradiation. A typical 260 W solar panel developed by LG Company 
was used for model evaluation using Newton-Raphson approach under 





Modeling MATLAB environment in order to analyze its behavior under actual 
Photovoltaic operating conditions. Comparison of our results with data taken from the 
manufacturer’s datasheet shows good agreement and confirms the validity of 
our model. Hence, the proposed approach can be an alternative to extract 
different parameters of any PV module to study and predict its performances. 
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1. INTRODUCTION 

Our environment is facing significant challenges in the present and future decades as energy 
consumption is expected to globally double during the first half of this century [1]. Due to increasing oil 
supplies constrains, the world is urged to search for alternative energy sources. Therefore, a better solution 
would be the investigation in Renewable Energy Sources (RES) [2-4] mainly, solar energy. Solar and other 
forms of renewable energies are seen as practical and clean solution to meet our planet’s growing 
environmental and energy challenges. In fact, the demand for solar energy is becoming more and more 
increasing and its consumption more and more generalized. The investigations on solar energy production, 
utilization and conversion are regularly increasing to cover more fields in the purpose of designing and 
implementing solar and electrical energy converters of higher performances and efficiency [5]. 

The most common principle widely known in solar energy conversion is the photovoltaic cell based 
on a p-n junction as a photodiode, i.e., it generates a voltage across its terminal when a light beam is incident 
on the device junction area [6]. Designers of electrical converters are interested in modeling PV cell/modules 
for studying the electrical converters including a PV system. This needs to know how to model the PV device 
feeding the converter. The PV modules show a nonlinear I-V characteristic with a number of parameters that 
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need be adjusted during experimental data analysis of practical modules. Mathematical modeling of PV 
cell/module is important to simulate the PV system and its components. 

Different methods for modeling photovoltaic (PV) cell/modules have been developed and explored 
earlier. Their main difference resides essentially in the electric equivalent circuit of the studied PV cell, their 
convergence, and complexity. Typically, Townsend [7, 8] model is one of the models developed for 
evaluating the parameters of the equivalent electric circuit of the cell. However, Rohan. S. Kolkarni and al 
proposed in their work another method for modeling the solar photovoltaic module using system 
identification [9]. However, Tarak Salmi and al have proposed a Matlab/Simulink model of PV cell based on 
mathematical equations using Simulink Bloks [10]. 

Zaid Hessein Ali et al [11] exposed a Simulation model of solar PV cell and hence PV panel using 
numerical approach consideraing only the effect of external parameters (temperature and solar irradiance) in 
western Iraq. The Matlab math modeling is implemented using climate and physical parameters with 
modeling equations [11, 12]. Newton-Raphson iteration numerical method has been applied to extract the 
value of current for every working voltage to find P-V curves under the effect of temperature and irradiation 
of Anbar province West Iraq [11, 13]. 

In this work, we demonstrate the behavior and functioning of a PV device by exploring its basic 
equations, modeling and simulating the 260 S1C-G2 PV Module using Newton-Raphson algorithm to solve 
the nonlinear I-V characteristic under Matlab environment. This is achieved by focusing on both the study of 
external parameters effects (temperature and irradiance) on the evolution of PV devices performances as well 
as those of the internal parameters (series resistance and parallel one). Comparison of the obtained data with 
that of manufacturer’s datasheet shows good agreement and confirms the validity of our model. 


2. PHOTOVOLTAIC GENERATOR 
2.1. Modeling the photovoltaic cell 

A solar module is defined as the individual piece of equipment that encompasses numerous solar 
cells connected in parallel or in series [14]. In order to simulate the real behavior of a photovoltaic cell, an 
equivalent electric circuit model is needed. Several models are proposed in literature such as the single diode 
model [15] which is derived from simplifications in the two diode model [16]. A single equivalent circuit 
model for a PV cell consists of a real diode in parallel with an ideal direct current source as presented 
in Figure 1. 

The two parameters used to model and characterize a PV cell are: the open circuit voltage (V,,) and 
the short circuit current (/,,). The V,, is the maximum voltage which a solar cell can provide at zero current. 
The /,,. is the maximum current which a solar cell can provide at zero voltage. 





Figure 1. Ideal PV cell model circuit diagram 


The output current from the PV cell can be found using the following equation 
I=- la (1) 


On the other hand, we have 
Isc = Iscrref (1 + a. (T; — Tres) (2) 
I,, at the reference temperature ( Iscrrer) is found on the datasheet, and a refers to the temperature 


coefficient of Ise in percent change per degree. The measurements are done under the standard irradiance of 
1000 W/m? and reference temperature of 25 °C. The photon generated current also varies with respect to the 
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irradiance G. The assumptions that /,. ~ Ip, is generally used since the series resistance (R, ) is low and the 
parallel resistance R,, is a very high. This leads us to the following expression 


lpn = r * Iscrref (1 + a.(T. - Trep)) 


G : Irradiation at temperature T, , Ger : Reference irradiance taken equals to standard irradiance of 1000 


W/m’, Tye Reference temperature equals to 25°C, And (Iq) is the shunted current through the intrinsic 
diode. Using Shockley’s diode model, the diode current can be written 


Iq = Ip. exp (sem) — 1} (4) 


A: Diode quality factor (it is between 1 and 5), K : Boltzmann’s constant (1.381.10°° J/K), T, : Junction 
temperature in Kelvin, q: Charge of the electron (1.602.10 “!? Coulomb), V : Voltage provided by solar cell. 
Ig is the diode saturation current, it depends also in temperature as it is expressed in relation (5) 


3 
IscT, Te \Y —qVgap /A-K 

I =— mn E) „exp | 212a (AK 5 

0 exp eee 1) Tref p Tai (5) 


AKT¢ Tc Tref 











I;cr,: Short-circuit current of the cell at temperature T. , Vjqy : Energy of the gap (1.16 eV for Silicon). 
Vocr, : Open circuit voltage at temperature T. , T;e¢ : Reference temperature equals to 25°C 
y : Similar to quality factor, it is taken equal to A. 

Several cells connected in series form the PV module and a collection of PV modules constitute a 
PV panel, and finally, the collection of a given number of panels forms the array. The practical current 
voltage I-V characteristic of a practical module requires the inclusion of the parameters (R, ) and (R,p, ) 
which makes the model more accurate as represented in Figure 2. 





Figure 2. Single diode PV cell model circuit diagram 


If we consider R,, = œ , we substitute the (2), (3), (4), (5) and (1) we get: 


3 
See IscTe T. \¥ —9-Vgap/AK V+ Rs) 
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(6) 

Taking into consideration the two parameters R,, and R, (6) becomes as follow: 
I Tt. \7 qep V+ RI V+Rg+1 
_ G A N scTe tel Y HAK + Rs. = te +Rs* 
I= Gref ` IscTe (1 + a. (T: Trer)) exp eTe) . (=) exp =a . exp {(a. a) 1} Ran 

Pe Te Tref 

(7) 


The series resistance (R,) accounts for the losses related to the connection of the cells in series, the 
resistance of the semiconductor material and that of metal grid whereas the parallel one (R;,;, ) accounts for 


Int J Pow Elec & Dri Syst Vol. 10, No. 4, Dec 2019 : 1914 — 1922 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 O 1917 





the losses associated to leakage current through a parallel resistance path to the device. It bcomes 
considerable for large number of parallel cells. These two parameters could be obtained experimentally from 
the I-V curve [16, 17]. This method is based on the fact that, the series resistance impacts significantly the I- 
V curve slope at nearby point (V,, ,0) and the parallel one at the point (0, I;..) , so we can write 


dv 


a (8) 
av 
R= = ia, (9) 
We have to derive (6) to get (R, ) expression 
Res Se o Lon o (10) 





ariy = 
V=Voc IscTc exp (t ocre) 
zap(4VocTc,)AKTe AKT¢ 
PUAKT¢ 


This model gives a better precision for many modules. The model expression results from simplifications 
obtained in the two diode model diagram presented in [16, 18, 19]. Cells connected in series increase the 
output voltage whereas cells connected in parallel provide higher currents values, so if our module is 
composed of N, cells then 


Iph(module) = Iph * Np and locmodule) =I * Np 


The current voltage characteristic of the PV module is a non linear equation; it should be solved 
using different methods. In this work, we have chosen the Newton-Raphson method for its quick 
convergence response as indicated in literature [16, 20]. 


3. .NEWTON-RAPHSON METHOD 

In this section, it is explained the results of research and at the same time is given the 
comprehensive discussion. Results can be presented in figures, graphs, tables and others that make the reader 
understand easily [2, 5]. The discussion can be made in several sub-chapters. 
Newton-Raphson algorithm consists in 


fxn) 
Xn+1 = Xn — Fay (11) 





f’ : is the derivative of the function f 
f(x) = 0: is the function to be solved,x, : is the present iteration, andx,,, : is the next iteration 
Then our equation to solve would be 








(1) = Ise -— 1 = Ip exp {(q * ss) -1}- ee =0 (12) 


The substitution of (12) in (11) we get 


z Patten) }- _Vtln Rs 
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4. RESULTS AND DISCUSSION 

In this section, the discussions of the results obtained in this work are divided into five parts. This is 
because the proposed model has been implemented and simulated under various operating conditions such as 
varied temperatures, varied irradiances, varying series resistances values and last varying the parallel 
resistances values. A comparison of obtained results with those taken from the manufacturer’s datasheet has 
been made. According earlier works presented in literature, the purpose of adjusting the mathematical I-V 
curve at the three perceptible points was successfully achieved. The electrical parameters of LG 260 S1C-G2 
PV module data are given in Table 1. 
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Table 1. LG 260 S1C- G2 PV module specifications under standard test conditions 
(1kW/m’, Am 1.5, T=25°C) 








Variable Value 
Maximum Power (W) Pmax 260 
Voltage at max power (V) Vupp 30,1 
Current at max power (A) Impr 8,64 
Open-circuit voltage (V) Voc 37,3 
Short-circuit current (A Isc 8,94 





4.1. Impact of irradiation 

Simulations are made at constant temperature of 25°C for various irradiations: (200, 400, 600, 800 
and 1000 W/m?), the Current Voltage and Power Voltage curves obtained are shown in Figure 4. From the 
results shown in Figure 4, it is obvious that the current module is directly proportional to irradiance as the 
photo current is directly proportional to irradiance and demonstrated in (3). It can be noticed that as 
irradiance value decreases, short-circuit current decreases proportionally. Decreasing irradiance also reduces 
the open circuit voltage (Voc), but following a logarithmic relationship leading relatively to a modest change 
at (Voc). Results illustrated in Figure 5 show that the module maximum power is also directly proportional to 
irradiance which makes the module as being more efficient when irradiance is important. 
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Figure 4. Irradiation level impact on current-voltage curves 
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Figure 5. Irradiation level impact on power-voltage curves 
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4.2. Impact of temperature 

Simulations are made at constant irradiance of 1000 W/m? for various temperatures: (0, 25, 50 and 
75°C), the Current Voltage and Power Voltage curves obtained are presented in Figure 6 and Figure 7. As 
can be seen in Figure 6, increasing module’s temperature decreases significantly the open-circuit voltage 
while the short- circuit current increases only slightly. In Figure 8, it can also be noticed that when 
the temperature of the module increases the power produced decreases, affecting negatively the 
module efficiency. 

The effect of temperature appears in the decrease of the open circuit voltage value. We conclude 
that maximum power generated by the PV generator depends strongly on the solar irradiation and 
temperature which means that the solar module generates a maximum power only for certain value of current 
and corresponding voltage. 
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Figure 6. Temperature impact on current-voltage curves 


Power Voltage curve for different temperatures 
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Figure 7. Temperature impact on power-voltage curves 


4.3. Impact of series resistance 

Figure 8 and Figure 9 show the simulation results of current-voltage and power-voltage curves 
obtained under the same Standard Test Conditions of temperature and irradiation for different series 
resistance values (R,). 

As seen in Figure 8, the series resistance of the module has a noticeable impact of the current- 
voltage curve slope in the region where the cell operates as a voltage generator. Integrating the series 
resistance effect in the PV equivalent circuit results in a voltage decrease at any given current. The increase 
in the series resistance value has a negative impact on the maximum power generated and hence, the 
efficiency decreases as illustrated in Figure 9. 
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Impact of series resistance value on the current voltage curve 
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Figure 8. Series resistance impact on current-voltage curves 
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Figure 9. Series Resistance Impact Power-Voltage Curves 


4.4. Impact of shunt resistance 

Simulations are made for Standard Test Conditions of temperature and irradiation (T=25°C and 
G=1000W/m7?). The equivalent electric circuit of a single diode model considered in this work is the one 
presented in Figure 2 considering R, = 8,9006. 1074 Ohm. 

Figure 11 and Figure 12 show the simulation results of current-voltage and power-voltage curves 
characteristic found under the same standard test conditions of temperature and irradiation and for different 
values of parallel resistance (R,, ). Figure 10 shows that the parallel resistance of the module has a great 
impact on the current-voltage curve slope in the region where the cell operates as a current generator. 

The increase in the parallel resistance value causes the increase in the value of the maximum power 
generated so the efficiency rises as demonstrated in the Figure 11. The shunt resistance serves as alternative 
paths for the free carriers produced by solar radiation. A higher shunt resistance means that a large amount of 
these carriers contribute to generate power whereas a lower shunt resistance indicates large losses, affecting 
mainly the slope of the I-V curve on the proximity of the short circuit region [21, 22]. The cell efficiency is 
affected because shunts reduce the Fill Factor (FF) and the Open Circuit Voltage (Voc). This effect becomes 
more prominant under low light conditions [23]. A low shunt resistance can lead to hot-spots in reverse cells, 
especially when the power dissipation occurs in a small area [24]. 
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Impact of Shunt resistance on current voltage curve 
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Figure 10. Parallel resistance impact on current- Figure 11. Parallel resistance impact on power- 
voltage curves voltage curves 


4.5. Comparison of results with the results taken from manufacturer’s datasheet and previous 
works modelling 

A comparison between calculated and experimental values of the considered PV module is made. 
The obtained results are presented in Figure 12. 

Results presented in Figure 12 shows consistencies between experimental and predicted results. It 
comes out evident that the calculated results are in good agreement with the experimental data provided by 
the manufacturer in datasheet. According to previous works presented in literature [10, 11, 15], the purpose 
of adjusting the mathematical I-V curve at the three perceptible points (Short Circuit Current, Open Circuit 
Voltage and Maximum Power) was successfully achieved. 
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Figure 12. Experimental and calculated I-V curves of the LG 260 S1C-G2 PV 
module under standard test conditions (25°C, 1000 W/m? and AM 1.5) 


5. CONCLUSION 

The represents work achieved a contribution to the prediction of photovoltaic device performances 
which constitute a very important step in the study of any PV system. The simplicity of a single diode model 
with the technique of adjusting the parameters and the improvements proposed in this paper can make this 
model an alternative tool for designers who are in need of a simple and effective model for simulating PV 
devices associated to power converters. So, the modeling simulation by Matlab is helful for PV systems 
designers due to the simplicity, effectiveness, accuracy and easy to use simulation and modeling method. 

The proposed model was used to study the impact of series and parallel resistances as well as 
temperature and irradiation on the development of photovoltaic device performances in terms of reliability 
and efficiency. Furthermore, simulation results show an excellent agreement with the experimental data taken 
directly from the manufacturer datasheet. Finaly our contribution could constitute a first step to build up a 
full solar PV power electronic conversion system in developing a grid connected application. 
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